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ABSTRACT 


A computer model of the launch and flyout of a shipboard surface 
to air missile toward an attacking aircraft or anti-ship missile has 
been developed for the Interactive Simulation of Engagements at Sea 
(ISEAS) program. The model, written in the C language, is based upon 
the MICE II Fortran program. The flyout model can be used to assess 
missile system performance against an air target. A complete discrip- 
tion of the missile model, with flow charts and a C program listing, 


is included. 
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I. INTRODUCTION 


ISEAS (Interactive Simulation of Engagements At Sea) is a computer 
program currently under developement at the Naval Postgraduate School 
to simulate interactions between a ship loaded with surface-to-air 
missile systems and an attacking aircraft or anti-ship missile. The 
portion of ISEAS described in this thesis is a translation of the 
Fortran program MICE II [Ref. 1:pp 4-28] for the flyout of a surface 
to air missile to the C language. Three different missile systems are 
simulated. The C program includes coordinate transformations, launcher 
bearing and elevation calculations, the proportional navigation equa- 
tion tor missile guidance, and the missile’s equations of motion. The 
program contains a description of the missile engagements with an air 
ΙΝ: ΠΡ ΡΡΟΡΙΕΞ Τα. Βα515 for the evaluation of surface-to-air 
missile launching opportunities and subsequent missile performance 
under realistic encounter conditions. This missile program provides 
the basic framework from which further simulations of increased comp- 


Ա enson isticutlom can be eosily implemented. 


Guided missiles are navigated along a path that leads to an inter- 


cept with 


II. PROPORTIONAL NAVIGATION 


the air target. Proportional Navigation is one type 


navigation in which the rate of change of the missile heading is 


proportional to rate of change of the line of sight between the missile 


and the target. See Figure 1 for the Proportional Navigation Guidance 


τι ἈΠΕ 
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O 
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missile heading 
(velocity vector) 


line of sight 


Figure 1. Proportional Navigation Guidance 
Trajectory 
Trajectory. The equation governing Proportional Navigation is 


where 


qe 


Qe 


8 = kg (1) 


time rate of change of missile heading 
(velocity vector) 


time rate of change of line of sight 


Proportional Navigation constant 


The missile rate of turn (8) is) G fixed proportion of the rate of 
türmsof theüne of sight (8), but the proportional constant itself (k) 
can be variable during the flight. This method is implemented in MICE 
τ" input base volwe κ. is multiplied by the ratio of the closing 
velocity to the missile velocity. Consequently, due to the changes in 
112 ο ο e value ot the Proportional Navigation constant 
is a variable during the missile flight. There are two special cases 
of Proportional Navigation, Pursuit and Constant Bearing or lead angle. 
The Pursuit Navigation course can be obtained from Proportional Naviga- 
tion by setting k = 1, and the Constant Bearing course can be obtained 
from Proportional Navigation by letting k be very large. 

A missile traveling a Proportional Navigation course is usually 
launched at a lead angle which corresponds to the Constant Bearing 
path. (Note that a pure Pursuit course is actually the unique case of 
a Proportional Navigation path with a zero lead angle and a navigation 
constant of one). Proportional Navigation paths are less curved (fewer 
g's on the missile) than the Pursuit course, but may be more curved 
than the Constant Bearing course. The advantage of the Proportional 
Navigation course is that it provides a practical method of causing an 
intercept with a manuvering target. Suitable values of the navigation 
constant are between three and six. |Ref. 2:p. 67) 

Figure 2 shows the two-dimensional geometry for a Proportional 
Γον οι ου course. The missile and target are a distance ro apart: 
The line of sight between the missile and the target is the angle 9 
relative to the global reference (X,Z) for all times of flight. The 


line of sight is at an angle €, relative to the missile velocity vector, 


and the line ot siqhthrsro on angle X, relative to the target velocity 
vector. The relative closing velocity between the missile and the 


target (r) can be written as [Ref. 2:p. 68] 


r= Ve cos, - VE cos OL ymn (2) 
where a cos c, - velocity components of the target along 


the line of sight 


V: cos X, = velocity components of the missile along 
the line of sight 


The rotation rate of the line of sight can be written in the 


form [Ref. 2:p. 68] 


rg = 2l Sin s V Sime. (3) 
where ան sin &, = velocity vector component of the target 


normal to the line of sight 


M sin X „= velocity vector component of the missile 
normal to the line of sight 


Ø = rate of change of the line of sight 
angle 


Equations (2) and (3) are the basic equations of motion for 


Proportional Navigation. 


A 


.- 


am, Via 
πα 


Figure 2. Proportional Navigation Geometry 
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ITIS 


A. GENERIC MISSILES 


THE MISSILE PROGRAM 


Simulations of three types of surface-to-air missiles; 


Medium Range, Short Range, are developed 


hor ras. 


Long Range, 


The general 


default missile parameters for each missile type are indicated below: 


Long Range 
wt (kg) 1066 
body dia (m) 0.55 
length (m) 7.98 
wing span (m) 8.91 
reference 

area ( sq. m ) 2778 
range (km) 121 
max alt (km) 20 
max flight 

time (sec) 200 
avg spd (mach) 2.5 
thrust (kN) 409 
initial 


velocity (m/sec) 670 


missile burn 
rate ( kg/sec ) 3.54 


drag coeff he 


11 


581 


@ . 55 


48 


20 


80 


409 


679 


Medium Range 


295 


Short Range 


251 


9.20 


5 ք 


1.00 


20 


20 


4 


200 


679 


67 


The user will have the option of changing any of the missile para- 
meters listed above, but after each missile firing run,  ISEAS will 


return to the original missile parameters. 


B. ENCOUNTER TERMINATION 
Ἱ.  Εεἰσπο πα Ծր 
Missile flight will be terminated if any of the following 
conditions are encountered: 
a. Missile range exceeds maximum range. 
b. Missile altitude exceeds maximum altitude. 
c. Missile time-of-flight exceeds maximum time of flight. 
2. Guidance Termination 
The ability to terminate missile guidance due to a target 


tracker break-lock is not incorporated. 


C. MISSILE POSITION IN FLCTCEHT 
1. Global Coordinate System 
The global cartesian coordinate system is used to define the 
target location. MICE II specified no reference point for its global 
coordinate, so for the ISEAS missile module the global coordinate is 
attached on sea surface with the Y vector to the north, the X vector to 
the east, and the Z vector perpendicular to X Tond yY Vector STOT CNE 
X,Y,Z coordinates (See Figure 3). 
2. Missile Stability Coordinate System 
The missile stability axis system is used to orient the missile 
in the global space, and its origin is at the missile’s center of 


gravity. The x axis is along the missile velocity vector (with 


12 





coordinate Systems 


Figure δ. 
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positive forward), the y axis is along the right wing, and the z axis 
is perpendiculor to the x and y axis (with positive downward). See 
Figure 5. 

The coordinates X,Y,Z locate the the missile’s center of gravity 
with respect to the global coordinates. The angles gamm and etom, 
shown in Figure 3, are the elevation and azimuth angles of the 
missile's velocity vector with respect to the horizontal plane of the 
global coordinates. Once the missile orientation is defined in the 
global coordinate system, transformation to the missile stability coor- 
dinate system can be performed using the following direction cosine 


matrix in the C language format. [Ref. 1:p. 4] 


qe 119112) τ[σ][1] T[8J[2] Ix 
ay f= | ոյ) ոյր TE22C23 Iy (4) 
ij Τ[21[6] τ[21[1] τ[21[2] Iz 


where 


Il cae the unit global coordinate vectors 


ΠΗ; - the unit missile coordinate vectors 


τ[6][0] 


cos(gamm)*cos(etam) 


ատն ը 


cos(gamm)*sin(etam) 


T[2][2] » sin(gamm) 


T[1][2] » sin(etom) 
1(1)(1) = cos(etam) 
T[1][2] = @ 
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IM 


-sin(gamm)*cos(etam) 


DR] 


-sin(gamm)*sin(etam) 


112131 


cos(gamm) 


For example, 


i cos(gamm)*cos(etam)*I, « cos(gamm)*sin(etam)*lys sin(gamm)*Ig 
is the unit vector of the missile stability coordinate in the direction 


of the missile velocity vector. 


D. LAUNCHER PARAMETERS 
A missile launch requires the following conditions: 


1. Line of sight range (from the missile to target) must be less than 
or equal to the maximum target lock-on range of seeker. 


2. Line of sight range must be less than the maximum missile launch 
range: 


3. Distance to the target after total missile flight time is less 
than or equal to the maximum missile range. 


ա ար 8 1-ի աար ԱՏԵ Gre Stationary and have no launcher bearing 
constroints. Since the ISEAS module's missile launcher is located 
ubpou5BdgESDEDpXuandsreolrstac Shipboard launch conditions are desired, the 
following limitations were added to the missile launcher requirements. 


4. Launcher bearings must not exceed 135 degrees port or starboard of 
snip s head! 


5. Launcher elevation must be less than or equal to 9Ø degrees. This 
will permit a vertical launch if desired. 


If any of the above requirements are not satisfied, missile launch- 
ing will not be permitted, and target tracking and intercept solution 
calculations will continue until all conditions are met. If launch 


conditions are satisfied, missile launch will be permitted. 


155 


The m_lchr function computes the launch angles (Azi and Elev) based 
on current ship and calculated target intercept positions. See Figure 
4 for the m lchr block diagram. Since the launcher is located on the 
ship, and the ship is in motion, launch angles will be relative to 
ship’s heading. The launch angles computation is interpreted from MICE 
II but with added launcher constraints as discussed above. The launch 


angle equations are [Ref. 1:p. 11] 


YTINT - YMSL 
Azi = Arctan[ ------------ ] (5) 
XTINT - XMSL 


ZTINT - ZMSL 
Elev = Arctan[----------------------------------------- TERG) 


sqrt(( XTINT - XMSL ) + ( YTINT - YMSL ) ) 


where 


II 


XMSL,YMSL, ZMSL the missile’s initial position 


XTINT,YTINT,ZTINT 


the target's position we ime ο 
intercept (provided by others). 


If the azimuth or elevation limits are exceeded, m_lchr will return 
to the main function in ISEAS for recomputation of intercept solutions 
or possibly for the user to maneuver the ship for better launch angles. 
If Azi and Elev are within the azimuth and elevation limits, the pro- 
gram will next check if the target is within range. If out of range, 
m lchr will return to the main fúnction and “ehe nadar wilteecont tame 


tracking the target. If the target is within range, m Kiriwi i em 
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Figure 4. M Lchr Flow Diagram 
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with the calculated azimuth and elevation values. This will be 
converted to angles in the global coordinate systemmond will οσο 


initial gamm and etam of the missile. 


E. MISSILE çi von RHA SE 

There are four functions that m flyout calls to compute the 
missile’s linear acceleration and angular velocities. These functions 
are the m atm(), m cna(), m thrust(), and m guide(). M flyout is 
discussed after these four functions. 

1. M Atm() 

The m atm function determines the atmospheric ΕΠΙ υπ 
temperatures as a function of missile altitude (ZMSL). MICE II employs 
a table look-up routine to determine the atmospheric temperature and 
density. The temperatures and air densities in  ISEAS are approximated 
by the following linear regression equations [Ref. 3:p. 299) based on 


the atmospheric data of Table 1 [Ref 1:p. 47]. See Figure 5 for m atm 


flow diagram. 


temp = 287.25 - 0.00624185*2MSL, for ZMSL < 11500 m 
(7) 
temp = 216.6, for ZMSL > 11500 m 
rho = 1.0177572 - 0.0000485359*Z72MSL (8) 
where 


temp - temperature of air at altitude ZMSL [ Kelvin ] 


rho = density of air at altitude ZMSL | kg/cu. m | 
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TABLE 1. ATMOSPHERIC DATA 


ALTITUDE TEMP CRATUPE PRESSURE OENSITY SPEED SOUND 


(FT) (D; K) (LB/SQ FT) (SLUG/CU FT) (FT/SEC) 
-2009.390 292.12 2273.72 0002519 1124.09 
0.09 283.16 2115,22 092377 1116.8 
2000.00 294.20 1967.69 002264 1105.7օ 
4030.00 230.24 1827.75 «002111 1100.99 
£000.08 276.29 1696.00 . 701947 1093.18 
8002, 00 Breese 1572.07 2001463 1045.32 
10000.00 263, 16 1455.60 . 031756 1077, 40 
12000.00 264-40 1345.23 «601648 1069. ԿՀ 
16090200 260.465 1243.6« 0001546 1061.39 
160C00.99 256.49 1147.49 00014449 1353.39 
14000.08 252.1 1357.47 «601355 1045.14 
20000. 00 240.57 373.27? 671767 1036.33 
22300. 00 244.52 494,59 «091154 1025.55 
24000.30 240.87 821.16 eGCULIS% 1329.39 
26000.00 235,71 752.71 0931029 1011.54 
24000.00 232.75 658,95 «002954 10063. «0 
30093, 00 228.41 629.66 «000891 994.85 
35300.30 218.93 493.3% «000754 973.1% 
0303.30 216.66 393.12 «000547 968.08 
45000.00 216456 309.45 «0326462 368.03 
50000.00 216.65 243.61 . 003368 368.05 
55000.00 216.60 191.50 «602257 964.03 
6300.00 215.66 151.03 «020226 368.04 
65303.00 216.66 118.933 «000174 963.08 
70905.99 216. 66 93,57 005160 365.048 
75000. 00 216. 66 73.78 «002117 968.08 
80Cc00.00 216.66 $8.13 «600057 963.08 
45000. 00 219.07 $5.83 «000064 973.4% 
90000. 00 223.542 35.29 . 009053 943.46 
95080.00 223.13 234.18 . 0000461 993.35 
102200.00 232. 65 ` 23.23 «000032 1903.19 
1656C2.20 237.19 15.54 . 001025 1212.31 
110055. 30 cum? 15.35 20020520 1322.52 
11902029520 246.26 12,10 «930916 1232.25 
220930 00 250.76 3.856 «020613 1051.44 
125090.00 255.28 8.93 . 000910 1050.42 
1360032, 89 259.80 5.57 . 070009 1C60.07 
140090. CO 258.82 4,06 «000035 1379. 13 
145633. 00 273.33 3.59 «211036 1927.34 
1$0000.90 277.86 3.05 . 001094 1035.27 
155006.909 282.35 2.55 2009093 1105.12 
160652.00 282.66 2.12 «000002 1105.74 
165022.00 282.66 1.77 «002002 1105.76 
176206.00 262.65 1.48 ». 000902 1105. Τι 
17 Stee. 60 282.66 1.23 «920001 1105.76 
130200.090 276.38 1.03 «000001 1033. 39 
185004. 20 259.65 . 85 . 009001 1979.94 
13109000 252.91 .70 ». 090601 1066.41 
{95900209 255.15 054 0930001 1252.55 
200900.09 249.45 «ο «009001 1238.76 
30953060490 C426 72 38 eC00221 1024.54 
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Figure 5. M_Atm Flow Diagram 
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Ze Chat.) 
ire Macnee tUneCE tone computes tne Slope of the curve that relates 
the lift on the missile to the angle of attack (the lift curve slope, 
cna) for all missiles [Ref. 2:p. 295]. MICE II used a table look-up 
for cna. In ISEAS, a different method is used. 
The lift curve slope is dependent upon the missile altitude. 


Ine sonic velocity at altitude is compúted using 


accel = sqrt( gamma * g * R * temp ) (9) 
where 
gamma = 1.4 
g = gravity const. [ kg-m/N-sq. sec ] 


R = 287 [ N-m/kg-Kelvin ] 


The missile mach number is determined by 


mach 


vmsl / accel (19) 
where 


vmsl 


missile velocity [ m/sec ] 


The lift curve slope coefficient is computed using the equation 


CHE b ΠΠ; 
where 


b = sart( mach * mach - 1 ) (12) 


See Figure 6 for m_cno flow diogrom. 
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Figure 6. M Cna Flow Diagram 
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5. M Guide() 

The m guide function computes the missile command angles, alpha 
(missile vertical angle of attack) and beta (missile side slip angle). 
See Figure 7 for the m guide flow diagram. 

The MICE II Proportional Navigation guidance equations are used 
in computing the angles, alpha and beta, which direct the missile's 
flight path. The solutions for alpha and beta are base upon the dif- 
ference in the target position between two consecutive time steps. For 
example, the closing velocity (closv) is determined by taking the 
difference of the range between the missile and the target from the 


previous time step, to the range of the current time step. Thus, 


closv = ( orng - rng ) / dt (15) 
where 
closv = closing velocity [m/sec] 
rng = new missile to target range [m] 
orng = old missile to target range [m] 
dt = time increment [sec] 


The direction cosines of the line of sight vector of the new time step 


in the global coordinate system were computed as follows: 


V[@] = ( XTAR - XMSL )/ rng in X direction 
v[1] = ( YTAR - YMSL )/ rng in Y direction (14) 
v[2] = ( ZTAR - ZMSL )/ rng in Z direction 


where 


XTAR, YTAR, ZTAR 


target position 


XMSL , YMSL , ZMSL 


missile position 
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Figure 7. M Guide Flow Diagram 
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The rate of change of the line of sight is computed using 


losdt[n] = ( Vial = vovin | ey dt (15) 


where 


losdt[n] rate of change of the line of sight 


for the nth direction 


V[n] » new direction cosine of line of sight 
oV[n] = old direction cosine of line of sight 
Mw. © 1.2 


Using the tranformation matrix presented in Section C, the rate of 
change of the line of sight is transformed from the global coordinate 


System to the missile stability coordinate system. 


mlil < mlil * t[i][j] * 1osdt[j] , m2 (16) 
where 
m[i] = rate of change of the line of sight 
along the ith direction 
t[i][j] = transformation matrix elements 
provided in Section C 
i cona Ø = x direction 
1 = y direction 
2 = z direction 
According to MICE II, alpha is given by [Ref. 1:p. 26] 
( k*closv*m[2] ) * W 
alpha = --------------------------------- ΓΙ»; 
cna * (5/2) * rho * vmsl * vmsl 
where 


m[2] = time rate of change of the direction cosine 
of the LOS vector in the z axis 


W - missile weight [kg] 


S = missile reference area [ sq. m | 
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Similarly, beta is οσο... 
(κ * closv * m[1] ) * W 
beta = -------------------------------- (18) 
cna * (S/2) * rho * vmsl * vmsl 


where 


mlL1] = time rate of change of the direction cosine 
of the LOS vector in the y axis 


5. M thrust® 
The m thrust function determines the missile thrust andemacsare 
weight as a function of time. MICE II used various numerical constants 
on its missile thrust and weight computations that were not thoroughly 
discussed. In ISEAS, missile thrust is taken as a constant value for 
each type of missile selected, and zero when the burn time is exceeded. 


The missile weight is computed using 


W = oW - ( burn rate * dt ) (19) 
where 
W = new missile weight [ kg ] 
oW - old missile weight [ kg ] 


burn rate - missile burn rate (depends on selected 
missile type) [ kg/sec ] 


dt = missile flight time increment 
See Figure 8 for the m_thrust flow diagram. 
6. M Flyout() 
Once the missile is fired by the user, them flyout 
function will check the following missile flight conditions at 


each time step: 
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Figure 8. M Thrust Flow Diagram 
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1. Missile altitude will be compared to the maximum missile altitude 
of 20,000 meters. 


2. Time of flight will be compared to the maximum missile flight 
time. 


3. Range between the missile and the target will be computed and 
compared to the maximum separation range allowed. 


If any of the above flight9ocheckscanevexcceededosthe mi flyout me eo 
will stop amd return to the main program. 

Another check is the minimun range between the missile and the 
target. If the range is 6.1 meters or less, this will be considered a 
hit. See Figure 9 for the m flyout flow diagram: 

Once alpha and beta are Computed by the m guide ην. ντ. 
conditions are generated and the linear acceleration equation for the 


vdotm is [Ref. 1:p. 6] 


[T - D - W * sin(gamm) - La * alpha - Lb * beta ] 


where 
T = missile thrust [N] 
D = missile drag at zero lift [N] 


La = cna τοσο ο 
missile lift due to angle alpha 


Lb = cna * q Տ * beto 
missile lift due to angle beta 
q = dynamic pressure = .5*rho*vmsl*vmsl [N/sq. m] 


See Figure 19 for the forces on the missile. 
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Figure 9. M_Flyout Flow Diagram 
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Calculation of the missile elevation angle turn rate (gamdm) 
due to the angle of attack, alpha, gravity and the lift component due 


to the net thrust is determine by using [Ref. 1:p. 6] 


[ (T - D) * alpha + La - W * g * cos(gamm) | 


Similarly the missile azimuth angle turn rate (etadm) due to 
the side slip angle, beta, is given by [Ref. 1:p. 6] 
| -(T - 0) * beta - ԼԵ * beta | 


etadm Ξ ------------------------------------------- (22) 
W * vmsl * cos(gomm) 





Figure 18. Forces On Missile 
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From the computed linear acceleration and angular velocities. 
an update of missile position, orientation and velocity is determined. 
In MICE II [Ref. 1:pp. 4-10], a slope averaging technique is used to 
numerically update the missile velocity and location. 


the missile is at initial position at time of fire. Thus, 


XMSL = XC vmsl = VC 
που πα ogamm = gamm 
ZMSL = ZC oetam = etam 


where 
XC,YC,ZC = values used in position corrections 
VC = values used in velocity corrections 
ogamm = old missile elevation angle 


Ὀανο SOUGemisslile azimuth angle 


At the first time step (dt), the missile is in flight and com- 


pe րթ ece linear Acceleration and angular velocities provide 


updates to missile position, orientation and velocity. [Ref. 1:p. 9] 
vmsl - VC + dt * vdotm (23) 
gamm = ogamm + dt * gamdm (24) 
etam = oetam + dt * etadm (25) 


XMSL = XC + .5 * dt * vmsl * cos(gamm) * cos(etam) (26) 
YMSL = YC + .5 * dt * vmsl * cos(gamm) * sin(etam) (27) 


ZMSL = ZC + .5 * dt * vmsl * sin(gamm) (28) 


3] 


The following are saved after the computation of equations (23) 
to (28): XC, YC, ZC, gamdm, etadm, vmsl, etam, gamm, vdotm, then are 
used as the old values, OXC, OYGC, OZC, ogamdm, oetadm, ovmsl, oetam, 
ogamm, ovdotm, for the next time step. 

As the missile continue on its flight path, at time greater 
than zero, gamdm and etadm values are generated and corrections are 


computed by the slope average technique. [Ref. 1:p. 9] 


GAMC = ogamm + .5 * dt * (ogamdm + gamdm) (29) 
ETAC = oetam + .5 * dt * (oetadm + etadm) (592) 
VC = oVC + .5 * dt * (ovdotm + vdotm) (51) 
XC = oXC + .5 * dt * (Cove + VC) * cos(GAMC) * cos(ETAC) (32) 
YC = oYC + .5 * dt * (oVC + VC) * cos(GAMC) * sin(ETAC) (33) 
ZC = oZC + .5 * dt * (oVC + VC) * sin(GAMC) (34) 
where 
GAMC = correction value for elevation angle 
ETAC = correction value for azimuth angle 


The missile position, orientation and velocity is updated by, 


[Ref. 1:p. 16] 


vmsl = oVC + dt * vdotn ΓΤ) 
gamm = ogamm + dt * gamdm (36) 
etam = oetam + dt * etadm Cay) 


XMSL = XC + .5 * dt * (VC + vmsl) * cos(gamm) * cos(etam) (38) 
YMSL = YC + .5 * dt * (VC + vmsl) * cos(gamm) * sin(etam) (39) 


ZMSL = ZC + .5 * dt * (VC + vmsl) * sin(gamm) (40) 
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The computed XC, YC, ZC, ETAC, GAMC,  vdotm, vmsl, etadm, and gamdm 
are saved and becomes the old values for the next time step. 

The output teem his functionewill consist of missile position 
(XMSL,YMSL,ZMSL), elevation angle (gamm), azimuth angle (etam), velo- 
ο ο (vist), weronteqw), sama missile thrust (1T). All of which the moin 


function can utilize for graphical presentations. 


F.o- KILL PROBABILITY 


This portion will conduct the endgame. It will be developed in a 


future study. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


The C programing language was used for this model because of its 
advertised user friendly dialogue techniques and its simple and effi- 
cient operations. Although the simplicity and generality favor the C 
programing language, since there is no standard definitions for the C 
language, compiling using Lattice C and Microsoft C compilers causes 
errors found from one unit but not from the other. It is therefore 
recommended to work with one selected compiler until the language is 
standardized. 

MICE II simulation logics was easily followed and understood. The 
translation to the C language for the thrust and weight calculations 
presented some problems because of the numerous undocumented constant 
variables used in the equations. The use of the table look-up method 
also created difficulties in translating, therefore, alternate equa- 
tions were used that provided basically the same results. All of the 
other major equations provided no additional difficulties in C language 
translation. 

The use of the missile model in decision making can assist the user 
by providing trade-offs and relationships and by evaluatiing and com- 
paring alternatives. But further work remains to be done. The missile 
model is currently working with only one guidance mode. Hopefully, ad- 
ditional guidance modes (command to line of sight guidance, beam rider, 


and track-via-missile) can be implemented. 
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Finally, during the early stages of this project, it is recommended 
that, (1) the student develope a working understanding on the mechanics 
of the equations of motion and guidance law, and (2) that the student 


complete a formal C language programing course. 
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APPENDIX A 


PROGRAM VARIABLES 


πας 2115 րա - ensino of computer program variables in 


alphabetical order. 


accel [m/sq.sec] 
alpha [rad] 
altmax [meters] 
azi [rad] 

beta [rad] 
ca 

el 


closv [m/sec] 
cna 


dt [sec] 
drag [N] 


elev [rad] 

etam [rad] 

etac 

etadm [rad/sec] 


g[kg-m/N-sq.sec] 
gamc 

gamm [rad] 

gamdm [rad/sec] 


ko 


lire [N] 
losdt [rad/sec] 


ոլ ] [m] 
mach 

mbrl [kg/sec] 
mbrm [kg/sec] 
mbrs [kg/sec] 
տուք 

m mr 

mselect 

IDs 


missile acceleration 

Missile attitude angle wrt stability axis 
maximum altitude 

relativer azimüth angle of missile launcher 


MES CNC çe Tip angle wre stability axis 


missile drag coeff 
missile lift coeff 
missile-target closing velocity 
missile lift curve slope coeff 


time increment 
missile drag 


relative elevation angle of missile launcher 
azimut angle or missile velocity vector 
correction value for azimuth angle 

rate of change of missile azimuth angle 


gravity 

correction value for elevation angle 
elevation angle of missile velocity vector 
rate of change of missile elevation angle 


base value of Proportional Navigation constant 


missile lift 
rate of change of LOS 


missile position in stability coord 
missile mach number | 

long range missile burn rate 

Πεν μμ ος πι πο ΟΙ rate 

short range missile burn rate 

long range missile (3) 

med range missile (2) 

missile select by user 

short range missile 
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oetadm 

oetam 

ogamdm 

ogamm 

orng [meters] 
ov[] [rad] 
ονο 

οχο 

ογο 

OZC 


q 


rho [kq/curmil 
rng [meters] 
rng1 [m] 

rng2 [m] 

rng3 [m] 
rtint [m] 


s [sq. meters] 


t [secs] 

t! [secs] 

t2 [secs] 

t3 [secs] 
temp [kelvin] 
thrust [N] 
thrust 1 oN] 
thrust2 [N] 
thrustš [N] 


v[] [radi 


νο 
vmsl [m/s] 


w [kg] 

wlr [kg] 
wmr [kg] 
wsr [kg] 


«Շ 
«ոտ1 [m] 
xtar [m] 
χει! 


yc 

ymsl [m] 
ytar [m] 
ytint [m] 


old rate of change azimuth angle 

old azimuth angle 

old rate of change elevation angle 

old elevation angle 

old range between missile and target 

old directicn cosine of LOS 

old correction value for missle velocity 
old correction value for xmsl 

old correction value for ymsl 

old correction value for zmst 


dynamic pressure 


atmospheric density 

range between missile and target 
max range of short range missile 
max range of med range missile 
max range of long range missile 
target intercept: ange 


surface area of missile 


missile flight time 

max flight time for short range missile 
max flight time for med range missile 
max flight time for long range missile 
atmospheric temperature 

missile thrust 

thrust of short range missile 

thrust of med range missile 

thrust of long range missile 


direction cosine of LOS 
correction value for vmsl 
missile velocity 


missile weight 

weight of long range missile 
weight of med range missile 
weight of short range missile 


correction value for xmsl 

missile position on x axis 

target post bere m << 

target intercept pest eircom tne x axis 


correction value for ymsl 

missile position on y axis 

target position ony Taxis 

target intercept position on the y axis 
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zc corcection value for zmsl 


zmsl [m] τ τ... Z axis (missile altitude) 
ztar [m] target position on z axis 
ztimtem target intercept position on the z axis 


53 


APPENDIX B 


PROGRAM LISTINGS 


The following is “Ες language programming ire em erme 


ISEAS missile simulaciomnm program 


Ն. 


m lchr() 


πι στο αἲι compute the louncher's azimuth and elevation 
melissa cm respect to the ship's heading. Input to this 
function will the launcher (the initial missile position) 
position and the calculated target intercept position. The 
output will be the launcher’s azimuth and elevation angles 
and must be converted to the global coordinates and becomes 
the initial gamm and etam for m flyout. “7 


1“ COORDINATES ARE RELATIVE TO SHIP'S HEAD */ 
ie time τπτ zLijt */“tgt intemeept positlon”/ 
fleGitax tare ycaor, ztLor- (“Lot rel position*/ 
τι ΠΠ} 251: /*loúncher or initial missile position*/ 
float elmts:=9ø.: /*max lchr elevation*/ 
ee mı 15 e << /*lchr port/stbd bearing Limits*/ 
frogat rng: range to target */ 
log i L s intercept range .*/ 
float rng1-20000,rng2-489009,rng521210090; /* missile ranges */ 
float elev, ozi; 
νοι χ.ν.Ζ.,ΧΥ' 


ieee CPU le LAUNCHER BEARING ANGLE eye / 
azi = atan((ytint-ymsl)/(xtint-xms1)); 
poe (182021 )/5.14159; 
rap: Ozi) > almts) 


{ 
ppodgtfo" ua BEARING ANGLE LIMITS EXCEEDED"); 
return; 
) 
2" ucNCoOMPERESMESSUNCHER ELEVATION ANGLE »*sswxJ/ 

x = xtint-xsml; 

y = ytint-ysml; 

ze ZL int-zsml 


Աա er yy) 

elev = atoan(z/xy); 

elev = (18@*elev)/3.14159; 
if (elev > elmts) 


{ 
printf("\n ELEVATION LIMITS EXCEEDED"); 
return; 

) 

e CHECK IF WITHIN MISSILE RANGE **<%**** / 
rtint = sqrt(x*x + y*y + z*z); 
if ( mselect == 1 ) 
{ 
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if ((ortint > rius ip 
( 
printf("\nTARCET NOTEWITHINT SEEKER RANGE"); 
printf("XnCONTINUE TRACKING"); 


return; 
) 

) 

if ( mselect == 2 ) 

( 

ifo ερ πι ὑπ. 
( 
printf("\nTARGET NOT WITHIN SEEKER RANGE"); 
orintf('"NncONTINUSE TRACKING"); 
return: 
) 

) 
if ( mselect == 3 ) 

( 

if (rtimn mance 
{ 
printf ("AnI ARGET NOT WITHIN SEEKER RANGE"); 
printf("\nCONTINUE TRACKING"); 
return; 
) 

) 


/**** PRINT LAUNCHER AZIMUTH AND ELEVATION ****/ 
lí az | 


orintf("\n LAUNCHER BEARING <2 2 (DEG το ορ 
else 


printf("Xn LAUNCHER BEARING 


%.2f DEG TO STBD" Gazi); 
printf("\n LAUNCHER ELEV = %.2f DEG",elev); 
printf("\n TARGET WITHIN RANGE"); 


orintf("\n LAUNCHER READY"): 
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finclude«stdio.h» 
finclude«math.h» 


/* This program was created to simulate the main function and 


αι ο τας cg ge ԱՅՐ ԱԱ ox voluate variable inputs and outputs 


ol the triyvoltrume e ron” 


Other functions are also called upon 


to provide the necessary values for the computations of the 


missile's equations of motion and eventually, 
Position in leva coordinates. 


the missile’s 


ա 


/* GLOBAL VARIABLES */ 


float xtarsz10090,ytar-10000,ztar-10000; /*current target coord*/ 
| loo men mi zasi. /“cürremt missile coord */ 

i om Yumsli-6 71 /* missile velocity - 2200 ft/sec */ 
float gamm,ogamm; /* elev angle of vmsl */ 

float gamdm, ogamdm; 

float etadm, cetadm; 

float etam,oetam; կոշ angle oi “vms 7 

float w: /* msl weight */ 

Iı eger. /* area of selected msl */ 

ilect 5152553: I ref area of short rng msl */ 

loat s2=4.92; /* ref area of med rng ms1 */ 
iibogtis5-9.78: /* ref area of long rng msl */ 

1 lom thrust: msi thirust 1 oN) rf 

i Koy eee cd ο, "msl drag.drag coeff */ 

eae roe., e msl lirt. lift coeff | 

float di: /* time increment delta-t */ 

float altmax-21000; /* msi maximum altitude in meters */ 
[uU EIC. /* msl flight time *7 

Please rag: 

float orng-20000; /* initial value from main function */ 
qplogtu5hostemp.9-d: 

Plog -alona, beta; 

roat cna, 

oat ko; /* input base value of proportional constant*/ 


int mselect; 


flog 
float 


միտ 125.|5|: 
ena; 


rrt shoot: 
int ons or si ansS2e,anss; 


froot 
float 
float 
float 
float 
froot 
TLlodt 
T LOGE 


YE XG WC σε 

(G po G FO yC OZC; 

vdotm,ovdotm; 

gamc,etac,.ogamc,oetac; 

ο ος mi Şi wir-1466; 

mors= 1 55 mbrm=0.558,mbrl=0.338; 

nd eee meets E ages 
thrust1-2000900,thrust2-400000, thrust35-5090090; 


43 


main( ) 


{ 
printf("\n Change target position? 1-yes/Ü-no "); 
scant("4d",. game. 
if ( ans == 1) 
{ 
printf( "VnEnter torget position im meters x 
թո ւոն ութը «եռ = ο. 
scanf("%f", &xtar); 
թո omre MEE 
scams "Xf yar) 
printf (“in ας ο Ὁ 


sconf("*f' &z tel 


) 
printf("\n Initial target position in meters. 
printf("\n xtar - 4f ytar = 41 Ztar = e e zi 
printf("\n Target heading is fixed ot 270 degs"); 
printf("\n at 18 meters decrement per time step"); 


printf("\n\n Enter missile initial position in meters"); 


printf (An xmms E .. 
scanf("4Zf", &xmsl); 
print? \m ymo 
scamf( gr gyms 
printf("\n zmn E 


scanf ( "Zf", &zms ih): 


printf("\n Enter missile attitude angle inlnadions E UR 
scant ("Sf"). &qeamm © 
printf( "Vn Enter missile dzimuthmangle an radions 
sconf("%f", &etom): 


II 
- 
- 

— 
-. 


printf("\n Desired time increment in secs = L 
Scan (EZ &dt 


printf("\n Enter proportional constant base value - Ho 
Οσο. δρ κο). 


printf('"wmEselectotypodor missmhew ο. 

printf("\nShort Range= 1 Med Range= 2 Long Range= 3 "); 
printf("\a missile select RS. 

scanf("4d", &mselect): 


printf("\n Change missile ref area? 1-yes/Ü-no  "); 
scan (“ran 
if ( ons) ==) 

( 
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1) 


if (mselect -- 


( 
printi ANM Enter: 51 = ye 
Տան ի ԱՏ RS]; 
) 
if (mselect -- 2) 
( 
premia Enter: s2 = n). 
ՏԵՑ ue Raso): 
) 
if (mselect == 3) 
{ 
pm Cm Enter: s3 = η): 
-επηί τν 255): 
} 
} 
printf("\n Change missile thrust? 1=yes/@=no աայ: 
Seamer. “sae &anse Ἱ' 
if ( ans2 == 1 ); 
{ 
if ( mselect == 1 ) 
( 
printi("Xn Enter: thrusti = e 
scan uz” &thrust 
) 
if ( mselect == 2 ) 
{ 
print © menter: tlmust? —.X); 
Sc If Si” wthrustay- 
) 
if ( mselect == 3 ) 
{ 
ΙΙ; nEErter:wtnpust»ses 7): 
Secom pes throust5», 
) 
) 
printf("\n Change missile burn rate? 1=yes/06=no "); 
Sean "sd". £anss5); 
if ( ans3 == 1 ) 
{ 
nim (mselect == 1) 
( 
bime m Enter: purn Թաետ 1.5. 71). 
scu mbrs)s 
) 
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if Լ mselece 2. 


( 
printf( "Vm Enter: puralraote 2 JN E 
scanf("%f", &mbrm); 
) 
if ( mselect -- 5 ) 
( 


կ 
E 
- 
W 


printf("\n Enter: Dorn ratë 3 
Scan ο emb 1 
) 


printf("ia When ready to shoot ο 
spacebar( ); 
t-o: 


wi Le toc uS 
( 
Orinitt (Anat =) oh et). 
m flyout(); 


«ԵԶՐԻ = «ԵԶԻ - 19: 
ορ πει πλω ίσο Sr xtar k 
ytar = ytar:; 
printeço ytor bp tor) 
ztor = ztar: 
printf” tup ն աի 
L. = İL. 
) 

print?é(' \neend οὐ ROC 

) 


spacebar( ) 

( 
printf("\n%s"," press spacebar"); 
while ( getch() !=' :)։ 

) 
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m flyout() 

/* This function computes the missile's linear accelerations 
musEGugulagr velocitics. Inputs of missile and target 
ΤΠ’ ΠΠ ΠΊΕ general characteristics are required. 
Complitat toms ot  missile-dttitüde and azimuth, velocity, 
Τι er ra αμ 111 600Ի6 the Outputs and are updated 
at each time increment iJ 

( 
(eros b 
tlogtaex yu 


x = xtar - xmsl: 
y = ytar - ymsl; 
Z = ztür >= zms1: 


(Դա: 38 96 3 HF range between missile and Carget KR E EEE E E EE / 


ON rt yy E z*z): 
printf("\n range = 4f meters",rng); 


[**"***xxxxx altitude check Of missile flight AKKM f 


if ( zmsl > altmax ) 
( 
printf("Xn max altitude exceeded "); 
return; 


) 


Yk CPA check of missile flight PEE IE HEDE DE IE IE DE DE GE JE IE IE IEE EEEE at at f 


e eOe Eo 2 --բչ 0061 km ör less */ 
return; /* warhead exploded- calc miss dist */ 


ար Արզն orno) 
( 
prine inimissile pass CRA"); 
“pap or nd. 
One meri = Zf m ",cpa); 
CERON 
return; 


) 


: Լն Ո function to calc temp and density wrt m) 
թ ՆԻ ա Սմ — return value of temp and density ******«** * / 


m atm(); 

örter (er anmaltitude of %f" 2msl); 
printf("\n temperature = 4f",temp); 
printf("\n density = %f",rho):; 
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[********* call m cna function to calc the lift curve slope E 
[********* coefficient ( cna ) with respect to altitude and * * / 
ΤΚ ΗΓ temperature «»/ 


mem ss 
orintf( "Xn emassa ՈՐՄ 


]******** call Ῥμπουίοπ ο σοἱ ο Ἡ Ὁ ο ουσ ιο ο ei 
9 return with thrust and msl weight E E o od o HE EE E re 9e Y Y e ve | 


dije ste 


[******* call function to calc msl guidance = return e 
/[******* alpha and beta — msi?commaudedgugles ο ο ο ο 2 
m guide(); 


[****xxxxwx calc drag at zero lift XOU e SEDE FE ICH EE t CO sss nn n 


printf mec η 

print (Anrhe e RO 
prime πι οἱ ΘΙ 
printf("Nm vmsl = %t" Umso 
drag » cd*(s/2)*rho*vmsl*vmsl; 
printf ("An droga aA E ὑπ 


[*******«*«** At time zero, missile is fire and initial +n ansa s 
[****"***»** misssile position is launcher position ԱԱ) 
if ( t es 9 ) 
( 
printf("Vn At cım one 
printf("\n xmsl= %f ymsl= %f zmsl- %f",xmsl,ymsl,zmsl); 
printf("\n etam= %f gamm= %f",etam, gamm); 


VC = vmsl; 

xG = «ms: 

yc s ymsl; /** values assigned to old values **/ 
zc - zmsl; /** for the next time step NO 
oetam - etam; 

ogamm - gamm; 

IE USE 

) 


[*******xxxwxxx Cale of msl linear accelaration XC X X ee o e x / 


GQ = οπαν ος πο vms πο  Ὁ μα 

b = σόπα”.οῦ 5 ρπο νησι υπο * beta; 

printf("\n a = @f b= gt ο. 
vdotm=((thrust-drag-(w*g*sin(gamm) )-(a*alpha+b*beta) )/w); 
printf("in vdotm = 4f", vdotm): 


48 


we Calc Of msl elavaction angle turn rate #363» iww / 


gamdm=( (thrust-drag)*alpha+a*w*g-(w*g*cos(gamm)))/(w*vmsl); 
printf("\n gamdm = %4f",gamdm); 


[ee = ~~ Colevoremsieozimuth angle turn rate *exx»xexJ 
etadm-(-(thrust-drag)*beta-b)/(w*vmsl*cos(gamm)); 


primet © mwcteddm Yam). 


a ie e emi ssile is in flight and if the time **/ 
το ΕΕΠ equals the time of increment, then first iteration**/ 


7 olimiss ii Tedpnosiciom orientation, and velocity is E 
ΤΗ computed. xD 
DER CN END 
( 


/* update of msl vel */ 
vmsl = ve + 9.5*dt*vdotm; 


/* new msl elevation angle */ 
gamm = ogamm + 2.5*dt*gamdm; 


/* new msl azimuth angle */ 
etam = oetam + %@.5*dt*etadm: 


/ supdace msl position base on flight correction */ 


xmsl = xc + 9.5*vmsl*dt*cos(gamm)*cos(etam); 
ymsl = yc + 9.5*vmsl*dt*cos(gamm)*sin(etam) ; 
zmsl = zc + 9.5*vmsl*dt*sin(gamm) ; 


/* change new value to old value for next iteration */ 


ovdotm = vdotm; 
ogamdm = gamdm; 
oetadm = etadm; 
ove = vmsl; 

OxC = xC; 

OyC = yC; 

OZC = ZC; 

oetam = etam; 
ogamm = gamm; 

) 
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[**** 
οσον, 
τ 
ՈՏԻՆ 


[ ***x* 


els 


/**** 


/[**** 


print meydan 


Since the time step is greater than the time increment **/ 
and greater than zero, the missile 15 ոմ c p I ww 
flight path and the missile position, ՓՈՐՆ ԱՆԿԱՆ s 
velocity are computed, and are corrected by the slope αν 
average technique աայ 


e 


gamc 
etac 
vc = 
xC = 
yc = 
ZC = 


vmsl 
gamm 
etam 
xms 1 
yms1 
zms1 


ov dot 
ogamd 


ogamm + .5 * dt * (ogamdm + gamdm); 


= oetam + .5 * dt * (oetadm + etadm); 


OVC 
Oxc 


ռս ո Ww ou 
N < X O O < 
o O 
+ O 
os 
3 3 


m = 
m = 


oetadm = 


oetam 
ogamm 
ove 
ος 
ογο 
οσα 


HHH HH 


HH HH OH OH 


+ dt * (ovdotm dom) 

dt * (ove + vc) * cos(gamc) * cos(etac); 

dt * Cove + ve) * cos(igame) ~ sin ecak 
է 


dt (ove + ve) * sin(gamc); 


+ + + 
տտտտ 
κ κα 


c + dt * vdotm; 

τ GOmam- 

+ dt * etadm; 

c+ .5 * dt * (vc + vmsl) * cos(gamm) * cos(etam); 
5 * dt * (vc + vmsl) * cos(gamm) * sin(etam): 
5 * dt * (vc + υπο) * sinfetom): 


vdotm; 
gamdm; 
etadm; 
Ετος, 
gamc ; 
VOS 

XC; 

yc; 

Ze 


Print missile position, orientation mand velocity mar 
as well as missile thrust and weight b 7 
printf("\nxmsl= %f m ymsl- Փո շոտ ԹԻ ո xmsl.vmsl,zmsT); 
primnti(” movmsi = “rm sects ums 


4f etam = Zf ",gamm,etam); 


printf("\n w = 4f kg πι το e 


retu 


) 


ՐՈ, 
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m atm() 
ΠΠ» function wilicalculgre the atmospheric density and 


temperature with respect to altitude. Equations used were 
derived by the linear regression method. Maximum missile 
altitude is 21000 meters, and the flyout subroutine will 
check if the missile exceeds maximum altitude. Therefore as 
long as the missile is at 21000 meters or below, atmospheric 
density and temperature will be computed. ԵՄ 
( 


[eee w density computation [ kg / cu. m | HHH HHH HHH HH HH / 


rhos( 1.0177572 - (06.00004585359 * zms1) ); 


[****** temperature computation [ kelvin ] "**x*****x**xxxxx**J/ 


if ( zmsl >= 11500 ) /** at this altitude, temperature **/ 
temp ~- 216.66: ա 5 constant Greene δο | **/ 
else 
temp - 287.25 - (0.006241865 * zms1); 


/** Return to flyout subroutine with values of temperature and 
density x 


91 


m cna() 

/* This function calculates che πε ου ο sov 
coefficeint with respect to altitude and temperature. CNA 
is used for the computation of alpha, beta, and the equations 
of motions κ) 


{ 
float accel: 
float mach: 
float b: 


(Դա calc missile aoceclordgi Tom -π 


accel = 20.05 * sqrt(temp); 


/** © Veale missile moch mumser աալ 


mach = vmsl / accel; 


b = sqrt((mach*mach) - 1); 


Որո Colee rit tmeurve stone ncocit xx f 
Շոց — 4 ο 


return: 
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m thrust() 

/* The m thrust function computes the missile weight during 
fiszgie-- Missile tnrust is G Constant value ond this value 
will depend on the type of missile selected. At burn out, 
missile thrust will be zero and flight time is exceeded. 


Output from this function will be missile weight and missile 
thrust. x 


/*** if user selects a short range missile ***/ 


ΜΙΤ -ΞΙ; 
{ 
Տ - si; 
if ( t == @ ) /** at missile time of fire κκ 
{ 
Wee WSr: 
thrust = thrust; 
return; 
) 
else 
( 
e I *xmissile flight time  **/ 
( 


w = w - (mbrs * dt); 
us = tnrustl1: 


return; 
) 
else 
( 
printf("\n missile burn time compeleted"); 
 .-.. 
πιεις Ր -.2: 
ΙΙ 
} 
} 
) 
/*** user selects a medium range missile ele | 
if Ար: - 155. -- ο ) 
( 
S = S2; 
.զ.-Ծ99Ք30Ծ:-ջ: missile time of fire κκ] 
{ 
ae ne 
E IRIS ea =" tir US to: 
return; 


9 


/ *** 


if ( mselect 


( 


else 


( 
if ( t « t2 Xx» missi leğpili ig mese ime 
( 
w = w - (mbrm * dt); 
thrust cans ee. 
return; 
) 
else 
( 
printf("\n missile burn time completed"); 
ο 
this ve vo 
) 
) 
if user select a long range missile x 
SN 
° = S5: 
if ( t == @ ) /*** missile time of fire ***/ 
{ 
w = wlr; 
tnrüst'scthrusts: 
return; 
) 
else 
( 
180ՆԵՀետ (2 == missile  tamemonme, tight = 3-7 
{ 
w = w = (mbr1 * dt); 
thrust = thrusts; 
return; 
) 
else 
( 
printf("\n missile burn time completed. 
w WwW: 
thrust = Ø; 
return; 
) 
) 
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m guide() 
/* This function computes the missile command angles, alpha 
(missile angle of attack) and beta (missile side slip angle), 


ο eoneropertciomplrevigation flight path. The command 
angles are with respect to the missile stability and are the 
output of this funecion: ui 

( 

arit cs 

ΠΠ I 

Todt mi Sil.: 


floot So 
floot Vosdt | oi]. 
ilogerclosy; 
LOOCK AN Zi 


x = xtar πο. 
y - ytar - ymsl; 
Z = Zeor.— zmsl: 
(moi sdPt( x*x + y*y + z*z ): 


ce XX omg ποτ: 


ն ՍԸ: ՀԱՐԹ nc of LOS in global coordinate ***/ 


առ = x/Zr ng 
ΙΙ 1 Բոց: 
Wei Zeno: 


/** calc the rate of change of LOS between missile and target **/ 


o We 52585: 12) 
( 
losdt[c] = ( v{c] - ov[c] )/ dt; 
"m une osdtl9d] s %f" c losdt[c]): 
) 


/*** matrix elements  ***/ 


աԱ = cosfgamm) * cos(etam); 
t[9][1] » cos(gamm) * sin(etam); 
L o pe ju Sam gamn); 
t[f1][9] » sin(etam); 


t[1][1] » -cos(etam); 

eliz 

ա ԱՆ Սարո): cos(etam): 
t[2][1] » -sin(gamm) * sin(etam); 


[ol IE scostaamm):; 
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/*** transform missile position from global coordinates to 
missile stability coordinates 
for ( i = ii <a 


ΚΝ 


{ 
m[i] = Ø; 
մազի... յ. 
m[i] > mba erene Tescile 
) 


/*** new direction cosine of LOS converted to old direction 


cosine of LOS for the next iteration. ας 
ον[6] = v[9]; 
ԽԱ νε. 
ov[2] = v[2]; 

/*** computed range between missile and target is now the 
old range for the next iteration. e 
rele a (ales 
printf("in old range = %f meters",orng); 

/*** calc the missile command angles, alpha and beta. թ» | 
alpha = ((ko*closv*m[2])*w) / (@.5*cna*s*rho*vmsl*vmsl1 ); 
beta - (ko*closv*m[1])*w / (9.5*cna*s*rho*vmsl*vmsl); 


printf (“n alpha = Ar Sna, 
printf( "Nn beta = Zf" Deta): 


return: 
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